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1. INTRODUCTION 


This report covers the work accomplished during the reporting period on the 
subject contract and describes plans for the next reporting period. The 
objectives during the current phase of the contract are divided into four phases. 
Phase 1 includes ground testing of the system and analysis of data from the 1981 
Severe Storms Test Flights. Phase 2 consists of preflight preparation and 
planning for the 1983 flight series. The flight test itself will be performed 
during Phase 3, and Phase 4 consists of post-flight analysis and operation of 
the system after that flight test. 

Section 2 describes the work performed during the present reporting period, 
and Section 3 provides plans for the next reporting period. 

2. WORK PERFORMED 

2.1 DATA ANALYSIS 

2.1.1 RANGE PERFORMANCE AGAINST niFFERENT TARGETS 

Figure 1 shows the range profile from five samples taken during Flight 10, 
around 1700 2. The lowest curve is taken from data collected upwind of Mt. 

Shasta at about 10,000 feet of altitude, in a clear atmosphere, where no signals 
were observed. It thus is a good representation of the noise level as a func- 
tion of range. The next curve was taken downwind of the mountain, and shows 
evidence of atmospheric returns. There is some question as to whether the 
data are valid at all ranges, or some ranges are contaminated by the others. 

The top three curves show the returns from Mt. Shasta at 8 to 10 kilo- 
meters. These demonstrate several features of the time history of the pulse. 
There is a "shelf" about 3 kilometers long, 9 dB below the peak of the pulse, 
and ringing after the pulse with a peak about 12 dB below the level of the 
pulse with a decay time of 33 microseconds and a period of almost 1 kilometer 
(6.7 microseconds) correspondi ng to a frequency of 150 kHz. The ringing does 
not decay to zero, but is still some 5 dB above the noise at the end of 200 
microseconds. 
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2.1.2 ANALYSIS 

OF PULSE 

SHAPE 






Det. 


Equi v . 

■ Rel . 

Rel . Energy 

Total Rel . 


Level 

Time 

Range 

Power 

per Gate 

Energy 


mV 

usee 

km 

dB 

dB 

dB 

Pulse Start 

142 

0 

0 

0 



End 

78 

2 

-0.3 

- 2.6 

0 

0 

Tail Start 

78 

2 

-0.3 

- 2.6 

- 2.6 


End 

0 

4 

-0.6 



-4.5 

Shelf Start 

0 

-78 

11.7 




End 

14 

0 

0 

-10.0 

- 9.0 

3.9 

Ringing Start 

6 

4 

-0.6 

-13.7 

-12.6 


Half 

3 

24 

-3.6 

' -16.8 

-15.6 - 


1/e 

2*2 

33 

-4.9 

-18.1 

-17.0 


Total* 






-4.0 

* Average over 







Peaks + Nulls 








Photographs of the pulse shape were taken by placing the Rockwell detector 
in front of the telescope with an aperture of a few millimeters and connecting 
the detector to an oscilloscope. The two channels of the scope were connected 
together to increase the gain, so that the pulse peak and the power before and 
after it were measurable. Several photos were taken with different vertical 
and horizontal scales and with different delays, to provide a time history of 
the entire pulse. These are shown in Figure 2. The center photo shows the 
entire time history. In this case, the pulse is barely visible in the third 
major box from the left. The lower two photos expand the regions before and 
after the pulse, and the top photo shows the detail of the pulse itself. 

The above table shows the results of analysis of the photographs of the 
pulse shape. The photos are considered to be divided into four parts. The 
first is the "shelf" which precedes the pulse. This occurs because the ampli- 
fiers are turned on prior to the optical pulse, in order to generate the 
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Figure 2. Pulse Shape Photos Taken at MSFC in January 1985 
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required gain. Unfortunately, this gain also amplifies the leakage tnrough 
the modulator. The shelf is described to a first approximation as being tri- 
angular, and the table shows the start and end of the shelf. 

The second region is the pulse itself, which is again described in 
trapezoidal approximation. The details of this region were taken from the 
top photograph. The third region is the tail of the pulse, before the 
ringing begins to be noticeable. 

Finally, the ringing region is considered. This region is caused by 
mechanical ringing of the modulator, and is characterized by sinusoidal modula- 
tion at a frequency of 150 kH 2 , which has a period of 6.6 microseconds, cor- 
responding to 990 meters of range. The envelope of the ringing is described 
by an exponential function, with the decay time calculated by measuring the 
time required to decay to half-height. The total energy is obtained by divid- 
ing the integral of the power envelope by two to account for averaging over 
many cycles. 

The immediate conclusions of this analysis are that the shelf can contri- 
bute to a spurious signal at a range before a hard target if the signal from 
that target is more than 9 dB above the noise. Furthermore, it will contribute 
a false velocity in this region if the hard target signal is more than 9 dB 
above the atmospheric signal from a range gate in the region. Finally, it 
will cause a similar effect in atmospheric signals if there are about eight 
or more range gates beyond the one being considered which have signals comparable 
to the current one. The total energy of 3.9 dB above the energy in the signal 
pulse means that in many cases, the signal will be corrupted by energy from 
this region. 

These conditions also hold for the ringing, except that the signals will 
be corrupted at ranges beyond the source of the signal rather than before it, 
and the magnitude of the effect will be smaller because of the lower powers 
and energies in this region. 

Application of these results to specific atmospheric measurements will be 
the subject of later work. 


Figure 3 shows the amplitude as a function of location for two forward 
and two aft scans. The flight path was from left to right, across the bottom 
of the plot. The ambient wind was from right to left at about 220 degrees 
relative to true north. It is evident that the atmosphere was clean and that 
no signals were detected upwind of the mountain. The strong mountain returns 
are visible about a quarter of the way from the right-hand edge of the 
pictures. Note that increasing darkness indicates increasing amplitude. The 
shelf before the pulse and the ringing afterward are clearly evident. Before 
the mountain signals are observed, there is an increase in the signal at near 
ranges, apparently due to increased backscatter in the air which has passed 
over and around the mountain. It appears that this occurs at progressi vely 
shorter ranges as the aircraft approached the mountain, and that the air became 
too clean to provide a signal, somewhere abreast of the mountain. The aft scans 
show a similar behavior. In addition, several scans are contaminated by strong 
ground returns at long ranges. 

Figure 4 shows the correspondi ng velocity plots. At the ranges where 
the intensities are enhanced, upwind (left) of the mountain, there are varia- 
tions in the velocity component which may be attributed to a wake caused by 
the mountain. When the mountain is visible in the scan, the velocity of the 
mountain is normally obtained at all ranges. This velocity has been examined 
in printouts of the data, and found to be within half a meter per second of 
zero. Lighter colors represent approaching (positive) velocities and darker 
ones represent receding velocities. The velocities were obtained from the 
VELOCITY data by changing the sign on the forward scans only, and adding the 
value of OFFSET. 


2.1.3 DECONVOLUTION OF RANGE DATA 

0 r e o* the major deficiencies of data collected during tre 1984 
flignl the existence of spurious data due to the failure of the 

electro-optic modulator. It is believed that a large number of the measured 
velocities are still valid, and that deconvolution techniques may be used to 
determine the data validity, thereby making most of this data usable for 
meteorological studies. 
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The problem may be summarized as follows: in normal operation, the lidar is 

a pulsed Doppler lidar, and returns processed in a particular range gate (a given 
time interval after the transmitter pulse) may be unambiguously associated with 
a range from the sensor determined by the given round-trip transit time. After 
the failure of the electro-optic modulator, a substantial amount of energy was 
transmitted both before and after the main pulse. Thus, in any range gate, the 
return signal is a sum of contributions from different ranges along the line 
of sight. In particular, at longer ranges where it is anticipated that there is 
no useful signal, contributions from the closer ranges, resulting from energy 
transmitted after the pulse, will often produce an apparent signal. The problem 
then is to determine how far out in range this signal is to be believed. 

The deconvolution techniques to be described below offer hope for resolving 
some of these questions. If full spectral data were available, it would be 
possible to perform a complete deconvolution and recover almost all of the lost 
information. However, in the system as it was used in 1984, only spectral 
parameters are available. Nevertheless, it is still possible to recover a 
considerable amount of useful information from careful evaluation of these param- 
eters as functions of range. The approach involves making assumption about the 
missing spectral information on the basis of available parameters. In this way, 
it is possible to establish varying degrees of confidence in the deconvolved data 
at different ranges. 

2.1.4 DERIVATION OF THE CONVOLUTION EXPRESSION 

The signal-to-noise ratio at any point in time after the transmitted pulse 
is an integral of contributions over all ranges from 0 to infinity. The average 
signal-to-noise ratio in a particular range gate is the average of the signal-to- 
noise ratio over that particular range gate. This can be represented- by the 
following equation, as shown in Figure 5. 
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Figure 5. Calculating the Contribution to a Range Gate 


10 



Raytheon 


where T is the duration of a range gate, 
n is the overall system efficiency, 

p( T j is the power transmitted, which is a function of the time at which 
it was transmitted, 
hv is the energy of a photon, 

B is the bandwidth, 

F(p) is the fraction of the transmitted power which is returned in such 
a way as to mix efficiently with the local oscillator, 
o is the range v and 

tg is the time associated with the center of the range gate. 

It will be noted that the two times and the range in the above equation are 

related as follows: 


The integral over range may be broken up into a sum of integrals over 
small segments. Later it will be convenient to make each of these segments 
correspond to a range increment equivalent to the time interval of the range 
gate. The resulting equation is shown here: 



where tg is the time associated with the center of the range gate, and 
R is the range associated with the center of the range increment. 

One may then conveniently define the average energy as follows: 

/ 

J T 

’ 2 



II 


and a weighted average si gnal -to-noi se ratio contribution from the range R as; 



2 2 


With this definition then, the measured signal-to-noi se ratio is related to the 
si gnal -to-noi se ratio contributions by the following equation, 

55r„ ■ E T 161 

n=0 

This equation relates the measured signal-to-noise ratio in a particular range 
gate to the contributions from all different ranges, weighted with the appropriate 
transmitter energy. It may be expressed in matrix form as follows: 

M = | C (7) 


In this case, Z is a column matrix in which each element is the signal-to-noise 
ratio contribution for a particular interval, starting with the closest and 
continuing to the furthest. Likewise, M is a column matrix whose elements are the 
measured signal-to-noise ratios in each of the range gates sequentially. The 
matrix £ consists of the elements E/E. Each of these elements is defined as shown 
below: 


m i 

* SNR m 

(tg * lAt) 

c j 

= SNR C 

R - J 

e i J 

II 

m| m| 

tt = (i - j) At 


2.1.5 CONVOLUTION 

The simple expression developed in the above section suggests the approach 
of inverting the matrix £, thereby solving for the SNR contributions in terms 
of the measured values; 
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C = £-1 M 


( 9 ) 


It can be shown that this is mathematically equivalent to the process of deconvolu- 
tion, which is often accomplished by taking the Fourier transform of the measured 
data and that of the convolving function, dividing the former by the latter, and 
taking the inverse Fourier transform. This provides an extremely straightforward 
method for extracting information where the intensity of a signal is a convolution 
along the line of sight of the desired function and the transmitter energy time 
history. All that is necessary is to calculate the matrix elements of E from the 
known transmitter time history. The resulting matrix may then be inverted and 
used to convert any measured data set to the desired output. 

In the case of the Doppler lidar, the convolution process applies to the 
signals obtained at each Doppler frequency. In the present system, complete 
spectral data are not available. The signal processor examines the spectrum 
which arises from the convolutions described above, and extracts the amplitude 
and Doppler frequency associated with the spectral peak. Thus, no information 
is available concerning the amplitude of frequencies other than that associated 
with the peak. One exception which will prove useful later on is the fact that 
it is known that all of these amplitudes must be less than or equal to the peak 
one. 

This fact may be applied in the following way. The matrix of measured 
data, M, for a specified Doppler frequency may be divided into two parts: the 
known part and the unknown part. The known part, M^, consists of elements whose 
value is equal to the measured amplitude in just those range gates in which the 
Doppler velocity rnder consideration is the one actually measured. The remain- 
ing elements are all zero. Conversely, the unknown part, My, consists of all 
zeros where the matrix has non-zero values, and vice versa. The remaining 
question is what values to use for these non-zero elements of M u . We will 
examine these in terms of their effect on the desired output. It is known that 
all of these non-zero elements must be non-negative since they represent 
measurements of energy in the receiver, and it is also known that the elements 
of E“lmay be either positive or negative. Thus the effect of an element of M u 
on an element of the output, C, may be to increase it or decrease it. The maximum 
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value of an element of C will be achieved when M u has large positive values in just 
those cases where they will be multiplied by positive elements of £, and zero 
•elsewhere. Of course, the maximum value which an element of M u may have is 
simply the measured signal in the corresponding range gate, since if it were 
any higher, it would become the new peak in the spectrum. This condition may be 
simulated by extracting from the positive elements, setting the negative 
elements to zero, in a new matrix The maximum value which an element of C 

may assume is then obtained by adding the value calculated by this matrix and 
the unknown measurements to calculate with the complete E*1 matrix and the known 
elements ; 


c T * r 1 M k + e; 1 m u (10) 

Likewise, the maximum value of an element of C is achieved when the unknown 
data points are at their maximum value in just those cases where they will be 
multiplied by negative elements of This may be achieved mathematically in an 

analagous fashion to Equation 11, with the following results: 

C N * I' 1 M k - E-l M u (11) 

Having performed the calculations indicated by Equations 9, 10 and 11, one 
obtains three values of the signal contribution for each velocity and each 
range gate. Comparisons among these values may be expected to yield a best 
estimate of the velocity at each range, along with at least a qualitative 
evaluation of its validity. 

2.1.6 TE' T RESULTS 

The procedures described above have been implemented on a few sets c r Jata 
from the 1984 Severe Storms Flight Test. Figure 2 shows the time history of 
the transmitter pulse. For mathematical modeling, the time history has been 
broken into four segments: the pre-pulse segment, the pulse, the tail, and 

the ringing. The energy in the time history in various times is given by 


pre-pul se 

E ■ _ 

in-0.9 | 

78 US + 

-78 us <t < 0 

E 

1 U 1 

pul se 

-§-■> 



t t - o 

tail 

-f ■ 10 

-0.26 


t t = T 

ringing 

o 

II 

|uj |uj 

- 1 . 26 e - 

(t t - t/29 us) 

o < t t 


H cos £ 

2tt 

(t t - t)J +1 



6.158 us 



and is computed by the FORTRAN function shown in Attachment 1. 

The deconvolution process is shown in Attachment 2. A portion of transformation 
E, and its inverse are shown in Figure 6. This figure shows in each case a 50 x 75 
segment of the 90 x 90 matrix, with the value of the matrix element encoded in the 
gray shade. It should be noted that in the display of the matrix itself, all values 
shown are positive while in the display of the inverse, both positive and negative 
values are shown. 

A set of data which included a measurement of Mt. Shasta was used to evaluate 
the deconvolution procedure. To account for difficulties in measuring the details 
of the ringing, a variable phase was included in the cosine term (Equation 12), and 
different values of phase were tried as shown in Figure 7. In each segment of the 
figure, the top curve illustrates the time history of the transmitter pulse which 
was assumed, the middle curve shows the measured signal as a function of range, 
and the bottom curve shows the result of the deconvolution. Note that an arbitrary 
time offset has been applied to the top curve. It is interesting to note that no 
significant variations seem to exist among the six figures shown. To show this 
in greater detail, three parameters of each figure were plotted as shown in Figure 8 
The variations among the different phases are small, and it appears that the best 
results are obtained at zero phase. In addition to the strong peak in amplitude 
indicating a return from the mountain at slightly less than 10 km, there is a 
smaller signal at close range, which may be atmospheric aerosol blown from the 
mountain into the beam, since the measurement was made slightly downwind of the 
mountain. After the mountain, a ringing signal still exists, since the decon- 
volution was apparently not accurate enough to remove this completely. 
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Figure 7. Mt. Shasta Measurements 
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The above example was particularly simple, since the velocity measurements were 
all nearly the same in a single deconvolution effort. To illustrate a more com- 
plicated example, typical of those expected in many of the interesting flow fields, 
a data set was chosen in which there was a distribution of velocities. In the set 
selected, most of the velocities did cluster around zero, -1, and +1 meters per 
second. In this case, Equations 9, 10, and 11 were all plotted in Figure 9. It 
is reasonably convincing from this figure that a velocity of zero meters per 
second existed in a region between about 6 and 9 km. Similar plots were constructed 
for the other two velocities. In addition, one other case needed to be considered: 
if there existed a velocity at which there was a significant spectral peak which 
remained just below the measured one at all ranges, it would be possible for this 
velocity never to appear in the deconvolution process. Thus, one case was run 
with % s These results are shown in Figure 10. In principle then, it is 
possible that the correct velocity for each range gate is never observed. However, 
one can make convincing continuity arguments that this situation is extremely 
unlikely. To better illustrate the partial deconvolution process in the case of 
multiple velocities along a given line of sight. Figure 11 shows several recon- 
structed spectra from the range sample shown in the previous two figures. Each 
figure shows four spectra. The top line, which is dotted, shows the peak of the 
envelope, indicating the highest peak which could possibly exist in the spectrum. 

The remaining three lines show from top to bottom the maximum expected and 
minimum signal at each velocity at the range written on the spectrum. At 7.5 km the 
velocity of 0 meters per second is very believable, due to the peak located there 
in both the minimum and expected spectrum. A similar, but slightly less convincing, 
situation exists at 1 meter per second at 6 km and again at -1 meter per second 
at 9 km. At 19.5 km, all the spectra are considerably lower than they were in the 
previous cases, and it is unlikely that it is possible to make a valid velocity 
estimate. 

2.2 MODULATOR REFURBISHMENT 

Upon receiving the two modulators back from the Marshall Space Flight Center, a 
number of tests were completed to recharacterize them. In the present Severe 
Storms system, the first modulator acts as an optical switch and the second as 
a "tail biter." During the last field test, the "tail biter" was electrically 
damaged, thus the need arose to reevaluate the modulator performance. 
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Figure 9- Three Convolutions at Zero Velocity 
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The extinction ratio and transmission of the modulators in tandem and also 
separately were measured. In addition to this, an infrared transmission scan 
was done. 

Upon completion of these experiments, the modulators were tested for optical 
ringing and the optical pulse height to optical ringing ratio was determined. 

The data from these two experiments indicated that two new crystals and modulator 
housings were needed. The crystals were ordered in April and a mid-November 
(previously mid-June) delivery is expected. The modulator housings were 
redesigned and fabricated to optical ringing. The design implemented has been 
successfully used on other Raytheon programs. Since the new crystals had such a 
long lead time, the original switching modulator was reworked (end faces repolished). 
This reworked crystal has been mounted in the new modulator housing and similar 
tests have been completed. The results of these tests will be summarized 
and forwarded to MSFC shortly. 

2.3 ISOLATOR REFURBISHMENT 

New indium antimonide was ordered and received in July. These wafers were 
polished by Cominco rather than the Research Division at Raytheon because they 
claimed they could do the work at minimal cost. The material was sliced from the 
original crystal grown in 1977. Examination of the received wafers indicated that 
Cominco did not have the skill in polishing as claimed. In addition, there is a 
surface mark on each wafer on the same face. (Crystals have an antimonide side and 
an indium side.) The surface mark is a sign of a material inclusion that propagates 
through the crystal and is therefore found in each slice. Inclusions can have a 
bad effect on the optical wave front as it passes through the crystal acting in 
effect as a lens since the index will be different than surrounding material. 

Tests are underway to examine the wafers at 10.6 pm. 

Old wafer stock, not yet polished, will be polished at Research as soon as 
blocking stock can be located. New stock will be ordered from Cominco in the 
raw cut state only. 
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2.4 IMPROVED CALIBRATION SETUP 


There are several problems with the existing internal calibration equipment. 
The signal-to-noise ratio (SNR) is very sensitive to the tolerance of the wheel 
placement due to the acute focusing of the present lens. In addition, the wheel- 
driving motor produces a non-optimum Doppler frequency, and the SNR exceeds the 
dynamic range of the processor. 

An existing 5.5X expander used in reverse would shrink the output beam sub- 
stantially to give a more narrow focal point. 

The wheel will continue to be placed at an angle of incidence of 45° due to 
the mass of accumulated data already processed at that angle. 

The motor, which drives the wheel, should be upgraded to raise the Doppler 
frequency (ffl) above 3 MHz. A 3.6K rpm motor would produce fp ■ 4.087 MHz. 

With a speed of 3.6K rpm, a beam size at the target of 4 mm would be 
required to produce a passband of 100 kHz. 

It has been determined that the optimum situation would produce a spot less 
than 1 mm on the wheel. Therefore, the length of time a particle is in the 
beam comes into effect to determine the passband. A particle on the disc (at 
the operating radius of 8.13 cm) will take 4.5 x 10~ 5 s to pass through the 
beam, which results in a 22-kHz passband. 

2.4.1 RESULTS AND DISCUSSION 

A diagram of the modification is depicted in Figure 12, showing the 
proposed 5.5X expander in reverse and dimensions of the system. 

The existing 5.5X expander design is recommended to be manufactured, and a 
3.6K rpm closed motor should be purchased. 

3. PLANS 

A new beam expander will be designed for the local oscillator. A cooperative 
effort with National Severe Storms Laboratory will begin on evaluating the 
velocity accuracy. 





At this point, all possible testing and fabrication have been completed. 
As soon as the new crystals arrive, they will be mounted in their respective 
modulator housings and individual and tandem testing will be finished. 
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